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ABSTRACT 

We describe the properties of the Lya forest in the column density range 
A^hi > 2 X 10^^ cm“^ based on 1056 lines in the wavelength range 4300 — 5100 A 
measured in extremely high S/N, R = 36,000 spectra of four quasars. The col¬ 
umn density distribution is well described by a —1.5 power law to 2 x 10^^ cm“^, 
below which limit confusion becomes too severe to measure a spectrum of indi¬ 
vidual clouds. The distribution of b values shows a well-dehned lower envelope 
with a cutoff at 5 = 20 km s~^ corresponding to a cloud temperature of 24,000 K. 

There is only a very small fraction (less than 1%) of narrow line clouds which 
cannot be identihed with metal-lines. From modeling the Lya absorption lines 
as complexes of clouds each with thermal broadening corresponding to be we hnd 
the b distribution can be understood if there is a mean of 3.25 clouds per ab¬ 
sorption line with a spread in velocity centroids characterized by a dispersion of 
10.75 km s“\ 

Subject headings: cosmology: early universe — cosmology: observations — galax¬ 
ies: intergalactic medium — galaxies: quasars: absorption lines — galaxies: 
quasars: individual (0014-1-813 = S5 0014-1-81) — galaxies: quasars: individ¬ 
ual (0302-003) — galaxies: quasars: individual (0636-1-680 = S4 0636-1-68) — 
galaxies: quasars: individual (0956-1-122) 


^Visiting Astronomer, W. M. Keck Observatory, jointly operated by the California Institute of Technology 
and the University of California. 
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1. Introduction 

The rich spectrum of Lyman-alpha absorption lines in quasars (generally known as the 
Lyman forest) was hrst noted by Lynds (1971) and studied in detail by Sargent et al. (1980). 
Because of its ubiquity and spatial correlation properties it is usual to assume that the forest 
is a distributed intergalactic population rather than being associated with large galaxies. In 
the hrst paper of the present series (Cowie et al. 1995) we analyzed the metallicity of the 
Lyman forest in clouds with column densities Nm ~ 3 x 10^"^ cm“^ and found that these 
objects have average metallicities of approximately 10“^ solar and consist of blends of much 
narrower components. In the present paper we turn to an analysis of the structure of lower 
column density forest clouds which, of necessity, must be based on studies of the Lyman 
absorption lines themselves since we do not currently have the sensitivity to detect metal 
lines in these objects. 

Prohle htting to high resolution spectra has been used to study the velocity width 
(Doppler parameter b) and the neutral hydrogen column densities N of the clouds. The 
distribution of Doppler parameters for the Lya forest was found to peak near 30 km s~^ with 
a rather large spread from 10 to 45 km s~^ (Carswell et al. 1984; Atwood et al. 1985; Carswell 
et al. 1987). Subsequently, on the basis of the very high resolution echelle observations which 
had become possible by the end of the 1980’s Pettini et al. (1990) claimed a tight correlation 
between b and N with many clouds having narrow prohles, leading to a median Doppler 
parameter of only 17 km s~^. This result created understandable excitement in theoretical 
circles since it seemed to question the standard picture of the intergalactic matter conhned 
in warm, extended, highly ionized gas clouds rather than small, cool, almost neutral entities. 
However, work by Carswell et al. (1991) using an identical observational setup did not conhrm 
the Pettini et al. conclusions. In a recent detailed analysis Rauch et al. (1993) showed that 
the prohle htting methodology is sensitive to noise and that this can produce the results 
seen by Pettini et al. 

The advent of the 10 m telescopes has resulted in an enormous improvement in the S/N 
which can be achieved in high resolution spectra of the faint quasars and in the present paper 
we discuss the distribution of low column density forest clouds and the b — N relationship 
based on snch data obtained with the HIRES spectrograph on the Keck telescope. We hnd 
that the b distribution is quite invariant as a function of column density and that there is a 
sharp cntoff in the nnmber of clonds below a valne of abont be = 20 km s~^, corresponding 
to temperature T = 24, 000 K for thermal broadening. We then argue that the broader 
absorption lines are blends of components, each with thermal broadening b^ and each broad 
line having an average of 3.25 components and a velocity dispersion of 10.75 km s“^. This 



structure agrees with that found from the metal lines in the absorption systems just above 
the present column density range. 


2. Data 

As part of a program to measure (D/H) ratios in high-;^ absorption systems, a number 
of bright (V < 18) z ^ 3 quasars with partial Lyman-limit edges (optical depth r < 3.5) 
have been observed with the HIRES spectrograph (Vogt et al. 1994) on the Keck 10 m 
telescope. A more complete description of the data and data reduction procedures is given 
in Songaila (1995). Each observation consists of a 3 — 6 hr exposure covering the wavelength 
range 3600 A —6000 A at a resolution R of 36,000. The wavelength coverage is complete at 
shorter wavelengths (A.$.5000A), but has small interorder gaps at longer wavelengths. Four 
of these spectra (toward the quasars 0014+813, 0302-003, 0636+680, and 0956+122) are 
well suited to a study of the low column-density (Vhi ~ 10^^ — 10^^ cm“^) Lya forest lines 
since they have good S/N (^50 per resolution element) and fully cover the wavelength range 
between the quasar Lya and Ly/? lines. The properties of the quasars are summarized in 
Table 1. 

For each of these four quasars a 600 A region of wavelength space was chosen which lies 
above the wavelength at which Ly/? would enter the spectrum but well below the quasar Lya 
emission wavelength to avoid proximity effects. The wavelength intervals are summarized in 
Table 1. The blaze function of the spectrograph in each order was removed by normalizing 
to hts to white dwarf standard stars obtained in the same conhguration before and after 
each quasar observation. The continuum level of the spectrum was then established by an 
iterative fit to each spectral order: all points above a critical value were htted by a third- 
order polynomial and renormalized; the critical value was then raised until the extrema in 
the noise in the Lyman forest region were matched to the extremal noise characteristics 
in the region above the rest-frame Lyman alpha emission taking into account the system 
response determined from the white dwarf observations and the overall spectral shape of the 
qnasar. The normalized spectra were then spliced together to form a total spectrnm. Each 
spectrnm was htted with an antomatic Voigt-prohle htting procednre which added lines at 
each minimnm in the spectrnm where the average absorption depth per resolution element, 
r, exceeded 0.05, nntil a satisfactory £t was achieved, in which the residnals were consistent 
with the expected noise levels. A sample region of a spectrnm with the hts and residnals is 
shown in Fig. 1. There were 1056 lines fonnd with iVni > 2 x 10^^ cm“^. The nnmber of 
detected lines for each qnasar is given in Table 1 together with the median b valnes in two 
column density ranges. The complete line catalogs and spectra are given in the Appendix. 
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For Q0014+813 we compared the detected clouds in detail with Rauch et a/.’s (1992) profile 
htting analysis of this quasar based on lower resolution (23 km s“^) and lower S/N data. 
The agreement was generally good and almost always within Rauch et a/.’s error estimates 
for both b and N. Almost all of the narrow lines in Rauch et a/.’s sample are found to be 
caused by scatter, which is consistent with their prediction. 

We have limited our primary analysis to the intermediate column density range 2 x 
10^^ cm“^ < A^hi < 3 X 10^^ cm“^ where uncertainties due to saturation are small. Measured 
column densities at Ahi > 3 x 10^"^ cm“^ are highly uncertain, but there are 66 clouds in 
the wavelength range with Ahi > 3 x 10^"^ cm“^. 

A plot of the b values of the htted lines vs. column densities is shown in Fig. 2 for A^hi 
in the range 3 x 10^^ cm“^ to 3 x 10^^ cm“^. The lower limit is chosen to avoid the optical 
depth selection effect against shallow lines, which is illustrated in Fig. 2 by the solid line, 
which corresponds to r ~ 0.05, and to minimize noise artifacts in the prohle htting. We 
have separated Q0636+680 (Fig. 2a), which has a relatively large number of narrow lines, 
from Q0014+813, Q0302-003, and Q0956+122, which have relatively few narrow lines. In 
Q0636+680 17 of the 266 lines in the 3 x 10^^ cm“^ < Ahi < 3 x 10^^ cm“^ range have 
6 < 14 km s~^. A detailed analysis of these narrow lines (Table 2) shows that 13 of these can 
straightforwardly be identihed as metal lines based on their correspondance to a strong Lya 
line and to other narrow lines in the spectrum. The remaining narrow lines are quite likely 
to be unidentihed metal lines. However, the most important point to draw from Q0636+680 
is that narrow lines can easily be retrieved from the spectrum if they are present, and that 
the absence of large numbers of narrow lines in the three other quasars is not a selection 
effect. The ease with which narrow lines can be picked out is also illustrated in Fig. 1. 

For Q0014+813, Q0302-003, and Q0956+122 (Fig. 2b) there are only 11 out of 670 
clouds with 3 x 10^^ cm“^ < Ahi < 3 x 10^^ cm“^ which have & < 14 km s“^. Eight of these 
are straightforwardly identified as metal lines (Table 2). Thus less than 1% of the Lya forest 
clouds have 6 < 14 km s“^. 

The b — N distribution of Fig. 2b has a fairly well dehned lower envelope which we have 
outlined with the dotted line 6 = 16 km s~^ +4 (Ahi/ 10^^ cm“^) km s~^. As we shall show in 
Fig. 3 the slope in this lower extremum is primarily an artifact of the much larger number of 
clouds at lower Ahi and the increased scatter in the b determinations, and the distribution 
is consistent with a sharp cutoff below 6 = 20 km s~^ which is invariant with column density. 


Turning to the distribution of column densities we show the raw distribution of counts 
per unit column density and 2 :, f d^dz \ quasars in Fig. 3. The counts are 



indistinguishable within the Poisson uncertainties, indicating that the properties of the forest 
are highly invariant and that the retrieval of the lines is reproducible from system to system. 


3. Analysis 

In analyzing the data of §2 it is necessary to deal with the problems of blending and 
selection since confusion becomes signihcant at the low column density end. Because there 
are a very large number of weak clouds the probability of cloud blending is high and can 
result in substantial incompleteness in the low column density counts and widening of the b 
value distribution. Very faint clouds may also be lost if their b values are too large and their 
central optical depths are too low compared to the noise levels. This latter selection effect, 
which is relatively unimportant at Ahi > 3 x 10^^ cm“^, is illustrated by the solid line in 
Fig. 2. 

In order to remove these effects we have performed a number of simulations to measure 
the cloud incompleteness and b distribution widening as a function of column density. In 
these simulations we assumed that the distribution of column densities was described by a 
—1.5 power law normalized to produce the observed number of lines, and that the distribution 
of b values was described by a Gaussian with b = 28 km s~^ and a{b) = 10 km s~^ truncated 
below b = 20 km s~^. We also assumed that the distributions were uncorrelated. These 
input parameters were chosen to produce output spectra which closely match the observed 
spectra. Artihcial spectra were generated and convolved with the instrumental spectral 
response with the appropriate Gaussian noise added. The hnding program was then run on 
these artihcial spectra to obtain lists of clouds. 

The completeness of the cloud recovery in the artihcial spectra is tabulated in Table 
3. Above 10^^ cm“^ essentially all clouds are recovered, together with a small number of 
blends of lower column density clouds, but below this column density the cloud recovery 
rate drops, falling to a completeness of only 25% in the interval Ahi = 2 x 10^^ cm“^ ^ 
4 X10^^ cm“^. It is important to emphasize that this ehect is not determined by the S/N since 
the incompleteness is almost wholly dominated by confusion: the recovery rate is essentially 
the same for an artihcial spectrum with inhnite S/N. 

The number of clouds per unit and unit Ahi in the observed spectra is shown in Fig. 3 
before and after correcting for the incompleteness. The corrected distribution in the range 
2 X 10^2 _ 3 

X 10^^ cm ^ is well ht by a power law of the form N ^ where N is the H l 
column density. The best ht to the power law index f3 is —1.46 with a 95% conhdence range 
of (-1.37 ^ -1.51). 



6 


Remarkably, the power law exponent of the spectrum htted over this limited column 
density agrees well with the power law £t determined over much larger iVn ranges (Tytler 
1987, Sargent et al. 1989). In Fig. 4 we have compared the present data with Petitjean 
et a/.’s (1993) summary of recent data at higher column densities and with a power law fit 
with (3 = —1.46. For consistency with their notation Fig. 4 shows the quantity 


f{N) = 


m 


AN 


( 1 ) 


where m is the number of clouds in the intervals dX and dN with X = 0.5[ (1 + zY — 1 ]. 
f{N) is roughly a factor of 3.8 less than ^ d^dz \ ~ present sample. 

We have tabulated f{N) in Table 3 for the present data. The two data sets agree extremely 
well in the overlap region. The present data set may also be used to check Petitjean et aids 
hnding that there is a significant deficit of clouds in the range 10^® — 10^^ cm“^ compared 
to the power law extrapolation from higher and lower column densities. The best power law 
fit to f{N) in the regime 2 x 10^^ — 3 x 10^^ cm“^, given by 


f{N) = 4.9 X 10^ 


( 2 ) 


with all quantities expressed in cgs units, would predict 193 clouds with Nm > 3 x 10^^ cm“^ 
compared to the 70 observed, implying that there is indeed a significant deficit between 
3 X 10^"^ cm“^ and 10^^ cm“^. Setting f3 = —1.51 only reduces the predicted number of clouds 
to 127. Thus there is indeed a significant steepening of the slope above Nm = 3 x 10^^ cm“^. 


The key issues for the b — N distribution is whether the rise of the lower envelope in 
b with N is real and whether there are any clouds with very large b or whether these are 
primarily random blends. We have compared the observed distribution of b values in the 
high (Ah I = 3 x 10^^ — 3 x 10^"^ cm“^) and low (Ahi = 3 x 10^^ — 3 x 10^^ cm“^) column 
density regimes in Fig. 5. The models give an acceptable fit to both distributions, showing 
that a single cutoff at & = 20 km s~^ can be adopted at all column densities. A better £t is 
given by adopting b = 18 km s“^ in the regime 3 x 10^^ — 3 x 10^^ cm“^ and b = 22 km s“^ 
in the 3 x 10^^ — 3 x 10^^ cm“^ range, so a weak increase in the cutoff b with N may be 
slightly favored but is not required. The drop in b with decreasing N may occur because 
weaker lines are more susceptible to having their Doppler parameters underestimated, due 
to the finite signal-to-noise of the data. 


4. Discussion 

Our first observational conclusion is that the number of Lya forest clouds rises as a 
smooth f3 = —1.5 power law in the range 2 x 10^^ — 3 x 10^^ cm“^ and shows no sign of 
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turnover at the low column density end. Confusion makes it essentially impossible to extend 
this result to much lower column densities. The second result is that the forest clouds have a 
well dehned lower envelope to the b values of 6 = 20 km s“^. It is possible for this limit to have 
a weak dependence on column density rising from b = 18 km s~^ at 3 x 10^^ — 3 x 10^^ cm“^ 
to 6 = 22 km s~^ at 3 x 10^^ — 3 x 10^"^ cm“^, but this is not required by the data. The 
distribution of b values can be adequately described by a Gaussian with b = 28 km s“^ and 
a a{b) = 10 km s~^ truncated below 6 = 20 km s~^ at all measured column densities. 

The mean intergalactic density, n, contributed by the clouds with these column densities 
may be obtained by integrating equation (2) whence 


Hn 


n = 


(l + x)=^ {l + 2qoz)2 


cNrr 


IN„ 


f{N)NdN 

m) 


( 3 ) 


Here F{N) is the mean ionization fraction at column density N. Assuming that F is not 
such a steep function of N that it weights the primary density contribution to the low column 
density end, we can integrate equation (3) to obtain 


fife 


1-5 X 10 ^ / Nmax 
< F > V3 X 10^^ cm"2 


( 4 ) 


where F is the appropriately weighted mean value and we have adopted Ho = 100 km s ^ 
Mpc“^ and go=0.5. The corresponding mean physical separation is 


L 


120 kpc 




2 X 10^2 cm-2 


0.46 


( 5 ) 


As Press & Rybicki (1993) have emphasized, equation (4) shows at once that the larger 
b values cannot be caused by the thermal broadening in the most simple photoionization 
model where the clouds are in thermal equilibrium and have no kinematic broadening. Thus 
from Donahue & Shull (1991) 

b = 28 km s"^ (6a) 

T = 4.9 X 10^ (6b) 

F = 3.4 X 10“® (6c) 

in terms of their ionization parameter U. However, from equation (4) we see that < F > 
must be greater than 1.1 x 10“^ in order not to exceed the baryon limit from nucleosynthesis 
of 0.01 (Walker et al. 1991) and hence U <0.3 implying b < 25 km s“\ 

The structure of the Civ lines at slightly higher column densities (Cowie et al. 1995) 
suggests an alternative picture. If we assume instead that each cloud is a blend of comparable 



column density components and that the observed b value is produced by this blend of 
components each with a thermally broadened be and a centroid velocity spread described 
by f{vc), then the value be must be identihed with the lower cutoff in b which will occur 
when an absorption line is dominated by a single component. This model is simplistic in 
assuming that all components have similar column densities and that there is a single thermal 
broadening for the clouds when in fact there are likely to be density and ionization parameter 
variations from component to component. However, the use of a single value for b is justihed 
by the extremely weak dependence of b on U in equation (6). Conversely, the extremely 
weak dependence of be on U also implies that U and F can be determined only to order of 
magnitude. Thus for be = 20 km s~^, U = 10“^ and F = 3.8 x 10“^, while for be = 18 km s~^, 
U = 10“^-^ and F = 1.7 x 10“^, and for be = 22 km s“^, U = 10“^-^ and F = 10“"^. 
However, and very reassuringly, this agrees with the U of 10“^ ^ 10“^'^ estimated by Cowie 
et al. (1995) based on the ionization balance in clouds just above Nm = 3 x 10^"^ cm“^. 
Furthermore, the hydrogen thermal broadening would predict a thermal b value of 6 km s~^ 
for carbon lines, which is very close to the value of 7 km s~^ found by Cowie et al. Finally, 
with F = 3 X 10“*^, then = 5 x 10“^ and very roughly 4% of the baryons inferred from 
nucleosynthesis he in these components. Thus this model appears to provide a self-consistent 
description. 

We can determine the average number of components and their velocity spread by 
comparing the predicted spread in b values with the observed distribution. If we assume 
that each absorption line consists of a number of components n, each thermally broadened 
with be = 20 km s~^ and described by the probability distribution 

P(n) ~ ^ (N)” with n > 1 (7) 

n\ 

and that the centroids of these velocity components are described by a Gaussian velocity 
distribution function characterized by o'(uc), then we can model the observed b value distri¬ 
bution function to determine N and Ve- We hnd that these quantities are extremely tightly 
constrained with a best fit of N=3.1, corresponding to a mean of 3.25 components per ab¬ 
sorption line, and a{ve) = 10.75 km s~^. Any smaller value of N results in too many single 
component clouds with b = be while larger N result in a Gaussian spread of b values which 
is also a poor £t to the observed distribution. The velocity spread is slightly smaller than 
the value of 18 km s~^ found by Cowie et al. (1995) for the spread in CiV components in 
Ah I > 3 X 10^^ cm“^ clouds and the mean number of Civ components (2.5) per absorption 
line systems is slightly smaller. This could represent a velocity dependent ionization effect 
or a difference between the column density samples, but most likely it is due to observa¬ 
tional uncertainty and reflects the difficulty of measuring the CIV dispersion. The presently 
determined dispersion should be considerably more accurate. 
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At larger velocity separations the components are no longer blended and we may directly 
investigate the two-point correlation function in the form 




^pairs ) 

(n) 

pairs \^/ 


-1 


( 8 ) 


where Npairs{v) is the number of observed pairs at separation v and N^^l^giv) is the number of 
expected pairs determined from the artihcial spectra with no correlation. This is illustrated 
in Fig. 6. For the forest sample with 7 x 10^^ cm“^ < Nm < 3 x 10^^ cm“^ we hnd no 
evidence of correlation at n > 150 km s~^ consistent with previous studies (Sargent et al. 
1980) but for 50 km s~^ < v < 150 km s~^ we do hnd a signihcant excess of 0.17 ± 0.045 - 
slightly smaller than, but consistent with, the value of 0.32 ± 0.08 found by Webb (1987). 
Intriguingly, if we divide the sample by column density the correlation increases with a value 
of 0.73 ± 0.13 for 4 x lO^^ < iVni < 3 x 10^^ cm-^ and 0.21 ± 0.07 for 7 x 10 ^^ < iVni < 
4 X 10^^ cm“^ suggesting that clouds are preferentially correlated with clouds of similar 
density, and that there might be a weakening of the correlation as we move to lower column 
densities. Crotts (1989), Chernomordik (1995), and Cristiani et al. (1995) have previously 
noted such an effect. Cristiani et al. found a correlation with a value of 0.89 ± 0.18 at 
Av < 100 km s~^ for Nm > 6 x 10^^ cm“^ clouds at slightly higher redshift in the spectrum 
of the quasar Q0055-269. For 2 x 10^^ cm“^ < iVni < 6 x 10^^ cm“^ they obtain 0.38 ±0.1. 
These values are roughly consistent with the present results, and together they may suggest 
that the correlation function does weaken at lower column densities. Suggestions of anti¬ 
correlation features in the two-point correlation at a few to several hundred km s“^ may 
also be noted in Fig. 6. While we cannot be sure of the nature of this phenomenon, such 
features appear real. Inspection of hgures 8 and 9 of Cristiani et al. (1995) shows that similar 
negative features are visible in the correlation function of their high column density sample, 
log IVhi > 14 (see also Meiksin & Bouchet 1995). Anticorrelation features are expected 
to occur when some Lya clouds are forming velocity caustics during turnaround (Kaiser 
1987; McGill 1990), but on smaller scales of a few hundred km s“^ they may be more likely 
to be due to gas dynamical effects like pressure fluctuations (Webb & Barcons 1991) or 
inhomogeneous ionization (Meiksin & Bouchet 1995). 


Finally, we may combine the present data with the properties inferred from the metal 
lines in the Nm = 3 x 10^"^ cm“^ Lya forest to estimate the cloud properties. Cowie 
et al. (1995) found from the ionization balance in the metals that U = 10“^ ^ 10“^'® for 
IVhi = 3 X 10^^ cm“^ clouds at this redshift, and, that the neutral hydrogen fraction was 
then 4.6 x 10“"^ (from Eq. 6(c) where U -2 = 77/10“^). By the dehnition of U the 

total density in the cloud is ritot = 6.3 x 10“^ J_2i/77_2 where J_ 2 i is the incident ionizing 
intensity just above the Lyman edge in units of 10“^^ ergs cm“^ s“^ Hz“^ sr“^. Assuming this 
ionization fraction holds over the relevant column density range, the line-of-sight cloud size. 
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H = 10 (-^Hi/10^^ cm“^). This is 10“® to 10“^ of the separation distance 

for TVhi = 10^"^ cm”^, suggesting that the hlling factor of the forest clouds is correspondingly 
small. The physical dimensions are similarly small compared to the transverse dimensions 
(Bechtold et al. 1994; Dinshaw et al. 1994), suggesting that the clouds are highly flattened. 
(Rauch & Haehnelt 1994 give a more indirect argument to this effect.) The internal density 

represents an overdensity of 1400 with respect to the average IGM, where VLjgm 

is the fraction of the closure density in baryons. 

The internal kinematic structure strongly favors models in which the underlying for¬ 
mation mechanism is gravitational or kinematic such as the minihalo model of Rees (1986) 
or the pancake model of Cen et al. (1994). The objects are extremely elongated, which 
would seem to favor the latter model. These models allow a natural understanding of the 
kinematic structure and the spread of velocities, since higher infall or shock velocities would 
result in too high a temperature for neutral hydrogen to be present. Indeed, the pancake 
model produces a natural explanation for the number of components, since we will see the 
pancake itself, the outgoing shocks, and possibly the surrounding infalling material. In its 
current version the pancake model does not correctly predict the extrapolation of the power 
law distribution of column densities below 10^^ cm“^ found in the present paper but this 
may be very much a function of the assumed details. 


APPENDIX 

Plots and line lists for spectral line hts for Q0014-I-813, Q0302-003, Q0636-I-680, and 
Q0956-I-122. Each plot shows the normalized spectrum with overlying prohle hts, plotted 
against vacuum heliocentric wavelength. Vertical tick marks indicate the positions of line 
features, with every tenth line numbered. Lines with Nm > 3 x 10^^ cm“^ and 6 < 14 km s~^ 
are marked by the hlled boxes. Residuals from the hts are plotted below each spectrum. A 
tabulation of the line identihcations by wavelength, with b value, and redshift and column 
density (under the assumption that these are Lya features) is given in Tables Al and A2. 


We would like to thank Jerry Ostriker, Jordi Miralda-Escude, and Martin Haehnelt for 
very helpful discussions. These observations would not have been possible without Steve 
Vogt’s HIRES spectrograph. 
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Fig. 1.— A sample 50 A region of the spectrum of Q0302-003 showing hts to identihed lines 
overlaid upon the normalized spectrum. Residuals to the hts are shown below, and it can be 
seen that there are no systematic residuals associated with the prohle hts. The wavelengths 
of each line are marked by ticks above the spectrum. There are two lines with 6 < 14 km s~^ 
in this region of spectrum where the ticks are overlaid by boxes. These lines correspond to 
a Civ doublet at z = 1.8924 (Table 2) and are easily distinguished from the broader forest 
lines. The number of lines htted for each quasar is summarized in Table 1. 


Fig. 2.— The distribution of components with respect to column density and b value for the 
quasar systems. Q0014+813, Q0302-003, and Q0956+122 have been plotted separately from 
Q0636+680, since the latter system has many narrow features known to correspond to metal 
lines. For each hgure we have plotted only clouds with 3 x 10^^ cm ^ < Ah I < 3 x 10^"^ cm ^ 
shown by the dashed lines. Very broad clouds with low column densities and high b values 
will not be detected. The solid line shows this selection for r = 0.05, above which clouds 
should be easily found. We have attempted to identify all cloud with 6 < 14 km s“^ as 
metal lines. Such narrow components that have been identihed with metal line systems are 
indicated with a surrounding box. The dotted line shows, for Fig. 2b, the apparent lower 
bound to the b values of Lya forest clouds. The boundary appears to be a real physical 
limit, and the plot of Q0636+680 shows that the htting program is quite capable of picking 
out narrow components where these are present. 


Fig. 3.— The raw and corrected distributions of counts per unit column density and red- 
shift in the column density regime 2 x 10^^ — 3 x 10^^ cm“^. Squares show the distribu¬ 
tion for Q0636-I-680, diamonds for Q0302-003, asterisks for Q0014-I-813, and triangles for 
Q0956-I-122. The dashed line shows the average raw counts for all four quasars. It can be 
seen that the raw counts are indistinguishable from system to system, within the Poisson 
noise, and that both the properties of the Lya forest and the capacity of the program to 
retrieve lines is highly invariant for all three quasars. The large open boxes show the column 
density distribution corrected for incompleteness with the best power law fit shown as a solid 
line. 
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Fig. 4.— Comparison of the cloud number distribution at low and high column densities. 
The data in the present survey, covering the column density range 2 x 10^^ — 3 x 10^^ cm“^, 
are shown with hlled squares, along with the best power-law £t (—1.46) over this region. This 
fit has been extrapolated to higher column densities, and the higher column density data 
compiled by Petitjean et al. (1993) is shown with diamonds overlaid with the associated error 
bars. The deviations below the line around 10^^ — 10^^ cm“^ which were noted by Petitjean 
et al. are confirmed by the present data (see text). However, as was originally emphasized 
by Tytler (1987), the overall fit to a single power-law is remarkably good over nearly ten 
orders of magnitude in column density. 


Fig. 5.— Comparison of the distribution of b values derived from the artificial spectra 
compared with the observed distributions. The lower left panel shows the input distribution 
for the artificial spectra averaged over eight realizations. In the right hand panels we show 
the recovered distributions (dashed lines) compared with the observed distributions (solid 
lines) in the high and low column density ranges. The observed distribution is the average 
for Q0014-I-813, Q0302-003, and Q0956-I-122 with known metal lines excluded. 


Fig. 6.— The correlation functions vs. velocity separation in the column density ranges 
4 X 10^^ — 3 X 10^^ cm“^ (upper graph) and 7 x 10^^ — 4 x 10^^ cm“^ (lower graph). The 
correlation functions have been computed in 50 km s~^ bins and the dotted line shows the 
la error computed from the number of artificial pairs. 



Table 1. Fitted Quasars 


Quasar 

m 

(1450 A) 

^em 

AA 

(A) 

No. lines 

N> 2 X 10^^ cm" 

No. lines 

-2iV> 3 X 10^^ cm-2 

Low N 

Median b 

(km s“^) 

High N 

Median b 

(km s“^) 

Q0014 + 813 

16.7 

3.384 

4500-5100 

254 

17 

30 

32 

Q0302 - 003 

17.8 

3.286 

4400-5000 

258 

12 

31 

32 

Q0636 + 680 

16.5 

3.174 

4300-4900 

302 

19 

26 

24 

Q0956 + 122 

17.8 

3.301 

4400-5000 

242 

18 

30 

29 


Note. — For each quasar, the number of lines htted and median b values are given for the 
indicated wavelength intervals in the observed frame. The number of lines htted by the program is 
given for systems with column densities in excess of Nm = 2 x 10^^ cm“^ and iVni = 3 x 10^"^ cm“^. 
Median b values are given for the column density ranges iVni = 3 x 10^^ — 3 x 10^^ cm“^ and 
A^hi = 3 X 10^3 _ 3 

X 10^^ cm Quasar magnitudes and emission-line redshifts are taken from 

Sargent et al. (1989). 
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Table 2. Narrow Line Identifications 


Quasar 

A 

(A) 

b 

(kms“^) 

ID 

Q0014 + 813 

4538.65 

9 

Cll (1334.5) 2.4009 


4582.39 

10 

— 


4597.20 

7 

— 

Q0302 - 003 

4478.01 

13 

Civ (1548.2) 1.8924 


4485.45 

11 

Civ (1550.8) 1.8924 


4832.47 

14 

Alll (1670.8) 1.8923 

Q0636 + 680 

4380.44 

7 

Sin (1260.4) 2.4754 


4390.19 

14 

— 


4417.85 

2 

Cll (1334.5) 2.3104 


4494.69 

10 

Civ (1548.2) 1.9032 


4537.95 

2 

Civ (1548.2) 1.9311 


4545.47 

7 

Civ (1550.8) 1.9310 


4545.81 

7 

Civ (1550.8) 1.9313 


4637.98 

3 

Cll (1334.5) 2.4754 


4642.63 

10 

— 


4646.95 

11 

Sim (1190.2) 2.9043 


4648.28 

3 

— 


4694.16 

10 

— 


4811.25 

2 

Mg II (2795.5) 0.7211 


4843.76 

9 

Civ (1548.2) 2.1287 


4848.16 

10 

Civ (1548.2) 2.1315 


4851.89 

8 

Civ (1550.8) 2.1287 


4856.22 

11 

Civ (1550.8) 2.1315 

Q0956 + 122 

4613.92 

6 

Si IV (1393.8) 2.3104 


4614.65 

8 

Si IV (1393.8) 2.3109 


4644.50 

8 

Si IV (1402.8) 2.3109 


4783.12 

11 

— 


4963.75 

8 

Sim (1206.5) 3.1141 
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Table 3. Column Density Distribution 


Logio N 

Range Incompleteness 

No. clouds 

Logio /(N) 

12.30-12.60 

0.25 

109 

-10.53 

12.60-12.90 

0.52 

188 

-10.91 

12.90-13.20 

0.82 

197 

-11.39 

13.20-13.51 

1.03 

173 

-11.85 

13.51-13.81 

1.11 

133 

-12.29 

13.81-14.11 

1.08 

88 

-12.76 

14.11-14.41 

1.01 

75 

-13.10 

> 14.41 

1.0 

77 

— 










b (KM/S) 





b (KM/S) 
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